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A class of autonomous quantum heat baths satisfying the eigenstate thermalization hypothesis 
(ETH) criteria is proposed. We show that such systems are expected to cause thermal relaxation 
of much smaller quantum systems coupled to one of the baths local observables. The process of 
thermalization is examined through residual fluctuations of local observables of the bath around their 
thermal values predicted by ETH. It is shown that such fluctuations perturb the small quantum 
system causing its decoherence to the thermal state. As an example, we investigate theoretically 
and numerically thermalization of a qubit coupled to a realistic ETH quantum heat bath. 
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I. INTRODUCTION 

Quantum thermodynamics has attracted much atten¬ 
tion in recent years. On the one hand, there has been 
much progress in understanding the origin of thermal¬ 
ization in closed quantum systems [1]. Thermalization 
is meant to be based purely on the phenomenology of 
quantum mechanics without any need for any external 
sources of heat to cause thermal relaxation of local ob¬ 
servables in a quantum system [2]. It was shown that 
it is the properties of eigenstates of an isolated quantum 
system whose local observables show thermal behavior 
[3114]. Later, a general canonical principle was intro¬ 
duced |S], stating that local observables thermalize for 
almost all pure states of a sufficiently big closed system 
under a global constrain. In this situation, the whole sys¬ 
tem with respect to local observables can be regarded to 
be in a microcanonical ensemble. On the other hand, in 
the studies of quantum thermal machines the emphasis 
is put on the properties of small systems which are cou¬ 
pled to some macroscopic sources of heat |6]. The heat 
baths are treated phenomenologically using the methods 
of open quantum systems. They are usually a collection 
of an infinite number of harmonic oscillators, which are 
presupposed to be in thermal equilibrium |7]. 

The smallest thermal machines are supposed to work 
at atomic scales and be of benefits in nanotechnology 
[Hj. To be able to do useful work they need to extract 
heat from heat baths. So far the role of heat baths has 
been played by macroscopic objects such as electromag¬ 
netic radiation or beams of atoms [5]. An interesting 
question arises whether it is possible to realize isolated 
quantum heat baths as sources of heat for such small 
machines. Combining together such small machines and 
microscopic heat baths it would be possible to create gen¬ 
uine quantum heat devices. A quantum heat bath must 
satisfy certain conditions to be able to thermalize a small 
quantum system coupled to it. Among them is the re¬ 
quirement that energies of the bath must be much denser 
than energies of the system and its density of states must 
be an increasing function of energy. Thermalization of 
the small system is then achieved by choosing a bath op¬ 


erator that couples to the system to be a random matrix 
m- The combined bath-system Hamiltonian is then it¬ 
self random with typically thermal behavior [5]. Indeed, 
it has been shown in Ref. m that a large class of isolated 
quantum systems in a pure state can serve as a heat bath 
for its subsystem. The crucial limitation of these works 
is that the bath and the system could not be separated; 
they were studied together as a closed quantum system. 
According to the second law of thermodynamics it is not 
possible to extract work from a single bath in a complete 
cycle; we need at least two baths to be able to build a 
thermal machine. The need for autonomous heat baths 
arises. In this paper, we present a class of finite real¬ 
istic quantum heat baths; those thermal properties are 
dictated solely by the properties of their eigenstates. We 
show that such systems, whose eigenstates satisfy the 
ETH criteria, cause thermal relaxation of smaller sys¬ 
tems weakly coupled to it. The ETH bath is thermal 
with some well-defined temperature regardless whether 
a small system is coupled to it or not. Therefore, the 
ETH bath is autonomous, the case which has not been 
addressed previously. This opens the possibility to use 
such baths to build thermal machines. ETH concerns the 
thermal properties of local observables of the bath, while 
the qubit is a separate system not covered by the ETH 
requirements. In this paper we show when and how the 
ETH bath thermalizes the qubit coupled to it. 

ETH states that under certain conditions almost all 
eigenstates of a closed quantum system with some Hamil¬ 
tonian Hb have the properties of a thermal state. Let 
|Efe) be an eigenstate of such system with eigenvalue E]^ 
and O its local operator. According to ETH if the di¬ 
agonal elements 0{Ek) = {Ek\0\Ek) are smooth func¬ 
tions of Ek, while the off-diagonal elements {Ek\0\Ei) 
are negligibly small, then for an initial state |(/>o) with 
energy Eb and small energy variance the expectation 
value {(j)\0\(j)) in the pure state |0) = exp{—itHB/h)\(j)o) 
relaxes to a steady thermal value and stays there for 
almost all times |3]. ETH predicts that the relaxed 
value tends to an appropriate microcanonical ensemble 

average (0|O|(/)) « ^ J2ki^k\0\Ek)■ Here we sum over 
JV energy eigenstates near energy Eb of the system. 
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Ek € [Eb — 5, Eb + S] with fixed 6 Eb- This yields 
that the whole isolated system with respect to a local 
observable can be regarded to be in the microcanonical 
ensemble with the density matrix 


Pm — \Ek){Ek\- (1) 

k 

For a finite quantum system, the values of the matrix ele¬ 
ments Oki = {Ek\0\Ei) are not ideally smooth functions 
of eigenvalues Ek but exhibit some residual fluctuations: 

Oki = 0{Ek)Ski + TZki, ( 2 ) 

where TZki are random with zero mean and small accord¬ 
ing to ETH. The semiclassical analysis shows that the 
values are characterized by a smooth function [T^ 

\TZki\^ :^SiEk-Ei)/2Trp(E), (3) 

where S{x) is the spectral function (see Ref. [12] for more 
details), p{x) is the density of states and E = {Ek+Ei)/2. 
The density of states grows exponentially with the num¬ 
ber of particles; therefore, the fluctuations are small if 
the number of particles is large. Recent developments 
have demonstrated that this behavior pertains to chaotic 
quantum systems lacking classical analogs [13] . The fluc¬ 
tuating behavior in Eq. ([^ manifests itself on the expec¬ 
tation value of the observable 

l^kl'^OiEk) + 

k k,l 

(4) 

The first term here is the expected relaxed thermal 
value, while the second term averages to zero at long 
times. However, the fluctuations {(j)\0'^\(f)) — {(j)\0\(j))'^ ~ 
Sfci are not negligible. They can be in¬ 

terpreted as residual thermal fluctuations [12] . 

II. ETH BATH 

Let an arbitrary small quantum system with Hamilto¬ 
nian Hs be coupled to a larger quantum system satisfy¬ 
ing ETH. For brevity, we call the latter ETH bath. The 
combined system Hamiltonian is given by 

H = Hs + Hi + Hb. (5) 

where Hi = gXs 0 O is the interaction between the sys¬ 
tem and the bath. Here g is an interaction strength; 
Xs is a system operator which couples to the bath op¬ 
erator O. Let |e;) be the eigenstates of the system and 
\Ek) be the eigenstates of the bath. The basis of the 


combined system is chosen to be \£ik) = \e-i) ® \Ek)- 
Using Eq. ([^, the combined Hamiltonian assumes the 
form of the sum of a “smooth” part with the elements 
{lk\H\Vk') = 5kk'[5w{^i+Ek)+gXivO{Ek)] and the “ir¬ 
regular” part with the elements gXuiTZkk'- Here, we de¬ 
noted Xw = (e/lAsje//). 

The smooth part is a block matrix which contains for a 
given index k a smaller matrix hiii(Ek) = 5ii'{^i + Ek) + 
gXu>0{Ek). The terms gXii'0{Ek) simply shift the en¬ 
ergies of the unperturbed system by an amount 5ei{Ek), 
which are some smooth functions of Ek- Eor an ETH 
bath the energy variance is small around mean energy 
Eb and we can assume that its eigenenergies Ek ~ Eb- 
Since the coupling g is small we may apply perturbation 
theory to estimate the energy shift 

5ei{EB) « gXuO{EB) + g^O^{EB) V (6) 

The irregular part, on the other hand, intertwines 
blocks with different indices fc in a chaotic manner due to 
randomness of TZw . This has a more profound effect on 
the dynamical evolution of the small quantum system. 
One can imagine that these random terms disturb the 
phase of the quantum state of the small system leading 
to its dephasing. Additionally, the off-diagonal nature 
of the irregular part allows energy exchange between the 
bath and the small system leading to dissipation. Both 
effects, dephasing and dissipation, lead to decoherence of 
the small system m- We will show that the ETH bath 
causes the system to decohere to a thermal state. 

HI. QUBIT COUPLED TO ETH BATH 

Consider a specific example of great interest: a qubit 
coupled to the above ETH quantum bath. The Hamil¬ 
tonian of the qubit with two available states |1) and |2) 
and corresponding energies ei and €2 reads 

Hs=ei|l)(l|+e2|2)(2|, (7) 

where A = £2 — ei > 0. We allow energy ex¬ 
change between the qubit and the bath considering 
for simplicity Xs = (t+ -I- (t_, where cr+ = |2)(1| 
and (T_ = 11) (21 are the raising and lowering oper¬ 
ators of the qubit. Following the formalism of open 
quantum systems [7], we move into the interaction 
picture Hi{t) = which 

yields ls(t) = and d{t) = 

As the next step, we need to 
calculate the correlation function of the bath operators 
C{t,s) = ii[0{t)0{s)pB\- As we discussed above, the 
bath satisfying ETH can be assumed to be in the micro- 
canonical ensemble pb = Pm, where the density matrix 
PttPS given in Eq. 0 . This allows us to evaluate the 
correlation function 
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FIG. 1. Elements of the density matrix of the qubit, which 
were obtained by solving Eq. (Ill and ( |12| | with the initial 
state pii = p 22 = pi 2 = P 21 = 0.5 and exp0A) = 4, 7 = 0.1, 
5A = 0.5, and A = 1. The qubit relaxes to the thermal state 
P 11 /P 22 = 4 and pi 2 = P 21 = 0. The process of thermaliza- 
tion (relaxation of diagonal elements) is twice faster than the 
process of decoherence (relaxation of off-diagonal elements). 


Cit, s)=C{t-s) = ^Y.T. 

k I 

( 8 ) 

The coupling is assumed to be small and the bath is 
large enough to satisfy ETH requirements. Under these 
conditions we use the Born-Markov approximation to 
study the dynamics [?]• In fact, it was shown in Ref. 
[14j that randomness of the coupling renders the second 
order approach valid. Under these assumptions we ob¬ 
tain the following Schrodinger picture master equation 


f^§iPs(t) = -i[Hs,ps{t)] 

poo 

-9^ dT[Xs,Xs{-T)ps{t)]C{T) 

Jo 

poo 

- 9 ^ dT[psit)Xs{-T),Xs]C{-T). (9) 

Jo 

To solve this equation we need to evaluate the integrals 
dTXs{—T)C{±T), which arise in the second and third 
lines of Eq. (1^. This boils down to taking the inte¬ 
grals I (A) — dTe=*=*^'^/^C'(T). By using the identity 

= TrS{e) ± where V denotes the principal 
value, we obtain 


/(A) = 5(A) ^ ^ p{Ek ± A)/p{Eu ± A/2) 

k 


+ 


N Ek-Ei±A- 


( 10 ) 


FIG. 2. Example of expectation values of the bath operators 
in the eigenenergy basis, (Efc|b°^b° = 0{Ek)5ki+TZki- At 
low energies the distribution of the diagonal elements 0{Ek) 
resembles smooth function in accordance with ETH, while 
at larger energies their behavior is more chaotic. Inset: off- 
diagonal elements Rk = \'Rkko\^ with fixed fco corresponding 
to Eko ~ 3.65Aeo. They are much smaller than the diagonal 
elements in accordance with ETH and their distribution re¬ 
sembles a smooth symmetric function around Ekg (shown in 
red) in accordance with the semiclassical estimate. 


Here, p(e) = ~ ^i) is the density of states of the 

bath and we used Eqs. ([^ and ([^ in deriving Eq. ( [Io| . 
The first term in /(A) is equal to 5(A)/2e^^^/^, where 
(3 = is the inverse temperature of the bath. 

We have assumed that A Eb, which allowed us to 
expand p{Ek ± A) « p{Ek){l ± /3A) « p{Ek)e^^^ and 
similarly for p{Ek ± A/2) « p{Ek)e^^^^'^. The second 
term in /(A) is significant if Ek = Ei when it is equal to 

Efc 0^iEk)/A « ±^O^iEB)/A. 

From Eq. (|^ we can now obtain two coupled equations 
for the diagonal elements of the density matrix describing 
the process of thermalization: 


= -g^5(A)e“^'^/^pii + 9^S{A)e^^/^p22, 

= -g^5(A)e^'^/^p22 + g^5(A)e“'^'^/^pii.(ll) 

The qubit relaxes to the expected thermal steady state 
limt_,.oo pii(f)/p 22 (t) = The specihc form of fluctu¬ 
ations \TZki\'^ used in the derivation of this result proves to 
be crucial to yield the correct thermal steady state. The 
remaining equations are for the off-diagonal elements de¬ 
scribing the process of decoherence: 


n 

n 


dpi2 

~W 

dp2i 


i{A + 26A)pi2 + i26Ap2i - "t[pi 2 - P 21 ), 
-i{A + 2 SA)p2i - i2dApi2 - 7(p2i - P 12 OI. 2 ) 


dt 
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FIG. 3. Inverse microcanonical temperature of the bath as 
a function of energy. It can be inferred from the entropy 
showing in the inset using the relation l3{E) = dS{E)/dE. 
Dashed line shows energy of the bath used in the calculations. 


where we denoted 6A = and 7 = 

(/^5(A) cosh(^A/2). It can be shown that the popula¬ 
tion relaxation occurs at the time scale 1 /( 27 ), while the 
loss of coherence occurs at the time scale I/ 7 . Thermal- 
ization of the qubit is thus twice faster than decoherence, 
which is a standard result in the context of atomic decay 
and quantum optics [7]. An example of this dynamics is 
shown in Fig. The off-diagonal terms can be shown to 
behave as ~ exp(it-\/A(A -I- 45A) — 7 ^//!) exp{—"ft/h). 
At small g we have 7 6 A A, therefore the 
energy levels of the qubit are shifted by the amount 
\/A(A 4 - 45A) - 72 - A « 25A = 20^{Eb)/A. This 
can be easily seen also from Eq. © , since in this case 
IA 12 I = 1 and £2 - Cl = A. 

We study now the case Xs = <^z, which does not al¬ 
low energy exchange between the bath and the qubit. 
Following the same procedure as above we arrive at the 
following master equations 

^ = —2g^5(0)pi2 -b iApi2, 

at 

^ = —2g^iS(0)p2i — iAp2i, 

at 

while dpxi/dt = dp 22 /dt = 0. This leads to dephas¬ 
ing of the off-diagonal elements \pi 2 {t)\ = \p 2 i{t)\ 
exp[—2(/^5(0)t/Ii] but not to thermalization, since pn 
and P 22 do not change. 

These results are in line with the analysis in the para¬ 
graph below Eq. Q: the first term in Eq. (§ causes 
energy level shift of the qubit, while the residual thermal 
fluctuations lead to the loss of coherence but not neces¬ 
sarily to thermalization (cf. also Ref. [E]). Therefore, 
energy exchange between the system and the bath is cru¬ 
cial to obtain thermalization of the system. 


The Markovian assumption might not be strictly satis¬ 
fied for a concrete system. In this case, for example, the 
upper limit of the time integral in Eq. ® cannot be ex¬ 
tended to infinity but only up to time t ^ . This renders 
more fine-grained dynamics especially at short times. As 
we are concerned with thermalization of a qubit, we are 
interested in the asymptotic state at t > 00 . It is thus 
expected that the above formalism is adequate for the 
present studies. 


IV. NUMERICAL EXPERIMENT 

To realize the proposed heat bath, we consider a two- 
band double-well potential filled with cold bosons. A 
complex interplay between the tunneling of bosons and 
their mutual interactions makes it possible to satisfy the 
ETH criteria and show thermalization m- We use this 
system as a heat bath for a qubit coupled to it. The 
Hamiltonian of the bath reads [13 [E] 

Hb = - - 1)+ J2 ^rK 

r^r\l r^l r,l 

+ c/°i (13) 

where h\. ( 6 / 1 ) are the annihilation (creation) oper¬ 
ators of a boson with mass niB in the well r = 
L,R on the energy level I = 0,1. For a given 

double well potential V (x) the corresponding single¬ 
particle wave functions (j)i.{x) can be found. The co¬ 
efficients in the above equation can be then evaluated 
[IB] as follows: The single-particle tunneling rate is 

j‘ = f dx(p’‘£{x) (^— 2 ^^^ + inter¬ 

action strength is = pB J dx\(j)l.{x)\‘^, single-particle 
energies are A/ = J dxcj^l*(x) (^-+ V{x)^ (j)l{x), 

and induced interaction between levels is = 

Pb J dx\(j)^{x)\‘^\(j)l{x)\‘^. Here, ps is the two-body inter¬ 
action coupling between bosons. 

The double well is created by splitting the har¬ 
monic potential vriB^QX^/2 by the focused laser 
lOfiwoexp(—a:^/ 2 (T^) with the width a = Q.l^JhJmBi^. 
We choose energy units eo = hujQ. For pB = 0.3eo we ob¬ 
tain jO = 0.26eo, = 0.34eo, U° = 0.14eo, = O.leo, 

= 0.06eo, Eo = 1.25eo, and Ei = 3.17eo. 

The Hamiltonian of the bath can be easily diagonalized 
numerically for N = 30 bosons and the eigenenergies \Ek) 
with the corresponding eigenvalues Ek can be extracted. 
The bath satisfies the ETH criteria at low energies as it is 
shown in Fig. We define the entropy of the bath at en¬ 
ergy E as S{E) = — Tr (pmlnpm) and the corresponding 
inverse temperature as its derivative /?(E) = dS{E)/dE. 
The results are presented in Fig. 

A two level qubit is represented via a single particle 
with mass ms trapped in a harmonic potential msuP'x^ (2 
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FIG. 4. Diagonal elements of the reduced density matrix 
of the qubit. It relaxes to the thermal state pwjp 2 i = 
exp(/3A). Solid lines correspond to the initial state of the 
qubit (1, l)^/\/2, dashed lines to (1,0)^ respectively. For our 
simulations we chose A = leo and we found j3 ~ This 

is in very good agreement with the inverse microcanonical 
temperatue of the bath; cf. Fig. Inset: corresponding 
energies of the bath. 



FIG. 5. The diagonal elements of the reduced density matrix 
of the qubit . The initial state of the qubit is (1, l)^/%/2 and 
the initial value of the bath energy Eb ~ 5Neo. ETH is not 
satisfied at this energy, cf. Fig. The qubit does not relax 
to the thermal state. Inset: corresponding energy of the bath. 


with two energy states ipnix) (n = 0 , 1 ) available to it. 
The corresponding single-particle Hamiltonian is 

1 

Hs = fiuj'^(n + l/2)aldn. (14) 

n—0 

We choose A = Huj = leg- The particle interacts with the 
bosons in the double well potential via contact interaction 
g. The Hamiltonian of the combined system is given in 


Eq. § , where the last term describes the interaction 
between the qubit and the bath: 

1 R, 

^1=9 Y. E (15) 

nn'll'—O r,r'—L 

where C"” ^ {x)4i^^i{x). The inter 

level transitions of the qubit between n ^ n' are allowed, 
such that the qubit and the bath may exchange energy. 
We expect the qubit to relax to a thermal state with the 
microcanonical temperature of the bath. 

We simulate the quantum dynamics by creating the 
Hamiltonian in the Fock basis of localized wavefunc- 
tions and propagate an initial state |'0(O)) by solving the 
Schddinger equation, \tp{t)) = exp{—iHt/h)\'tlj{0)). The 
initial state of the entire system is a Fock state of the bath 
|n° ,n^,n];^,n)j) times an initial state of the qubit. We 
choose two initial states for the qubit: ( 1 , 0 )^ correspond¬ 
ing to the particle occupying initially the lowest state 
of the trapping potential and (1,1)^ /v^ corresponding 
to the superposition of the lowest and the first excited 
states. The initial Fock state of the bath is chosen such 
that the energy of the bath {iIj{0)\Hb\'4’{^)) ~ 3.65fVeo 
satisfies the ETH. The energy of the bath 
changes in time slightly from this value, since the qubit 
is coupled weakly to the bath and its Hilbert space is 
much smaller than that of the bath. At this energy the 
inverse temperature of the bath can be found from Fig. 

/3 « O.Scq As expected, the reduced density matrix 
of the qubit p{t) = TrB\'tp{t)){'ip{t)\ relaxes to the thermal 
state limt_,,oo pii(t)/p 22 (i) = exp(/3A) as shown in Fig. 

On the contrary, thermalization is not observed for 
higher energies of the bath where ETH is not satisfied as 
it is evident from Fig. Therefore, ETH is crucial to 
obtain thermalization of a small quantum system. 

V. CONCLUSIONS 

We have demonstrated that isolated quantum systems 
satisfying the criteria of ETH can serve as finite au¬ 
tonomous heat baths for smaller quantum systems. As 
an example, we studied theoretically and numerically the 
thermal relaxation of a qubit weakly coupled to a realis¬ 
tic ETH heat bath. We believe ETH heat baths will be of 
benefit in building thermal machines working genuinely 
on the microscopic level, where not only the engine but 
also heat bath is treated quantum mechanically. 

The qubit is believed to be the smallest quantum en¬ 
gine. It is an interesting open question of how large the 
Hilbert space of the thermalized system can get. 

The author acknowledges support from the Marsden 
Fund (Project No. MAU1205), administrated by the 
Royal Society of New Zealand. 
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